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We develop a proper nonempirical spin-density formalism for the van der Waals density functional 
(vdW-DF) method. We show that this generalization, termed svdW-DF, is firmly rooted in the 
single-particle nature of exchange and we test it on a range of spin systems. We investigate in detail 
the role of spin in the nonlocal-correlation driven adsorption of H 2 and CO 2 in the linear magnets Mn- 
MOF74, Fe-MOF74, Co-MOF74, and Ni-MOF74. In all cases, we find that spin plays a significant 
role during the adsorption process despite the general weakness of the molecular-magnetic responses. 

The case of CO 2 adsorption in Ni-MOF74 is particularly interesting, as the inclusion of spin effects 
results in an increased attraction, opposite to what the diamagnetic nature of CO 2 would suggest. 

We explain this counter-intuitive result, tracking the behavior to a coincidental hybridization of the 
O p states with the Ni d states in the down-spin channel. More generally, by providing insight on 
nonlocal correlation in concert with spin effects, our nonempirical svdW-DF method opens the door 
for a deeper understanding of weak nonlocal magnetic interactions. 

PACS numbers: 71.15.Mb, 31.15.ej, 81.05.Rm 


The modular building-block nature of metal organic 
frameworks (MOFs) and their extraordinary affinity for 
adsorption of small molecules make these nano-porous 
materials ideal for technologically important applica¬ 
tions. MOFs are used, for example, for gas storage 
and sequestration catalysis IS polymerization 

ISIS], luminescence [TOl [TT] , non-linear optics [12] , mag¬ 
netic networks [13] , targeted drug delivery [14] , multifer- 
roics nsHni, and sensing [T8142T] . The design of novel 
MOFs with improved properties requires insight into the 
molecule/MOF interaction. The large unit cells and 
periodic nature of MOFs make density functional the¬ 
ory (DFT) the prospective tool for a theory exploration. 
However, both the adsorbate molecule and the MOF’s 
metal centers can carry spin, giving rise to complex mag¬ 
netic interactions and a molecular-spin response. It is 
thus crucial that DFT can reliably capture van der Waals 
(vdW) forces—which govern adsorption in MOFs—in 
concert with spin effects. 

Concerning the former, the last decade witnessed the 
development of DFT descriptions for these forces [22]. 
Here, the vdW-DF versions [231126] stand out by be¬ 
ing nonempirical exchange-correlation functionals that 
are systematic and truly nonlocal extensions beyond 
LDA [27] and GGA [28] in the electron-gas tradition 
[22] [29] [30] . Subsequent developments include variants 
which differ by their choice of the semi-local exchange 
[3TII35] and related nonlocal correlation functionals that 
rely on optimizing parameters [36H38] . The vdW-DF 
method and relatives have been successfully applied to 
numerous materials in general isaEniisn], and to small- 
molecule adsorption in MOFs in particular [l|5l|40H46]. 

Goncerning the spin effects, however, a systematic de¬ 
scription within the vdW-DF framework is still missing. 
Such effects can play important roles not only in MOFs, 


but in many systems, as Hand’s rules reflect a prefer¬ 
ence for spin-polarized ground states. For example, spin 
and vdW effects are essential in organic spintronics m, 
dimer binding in excited states [48], overlayer forma¬ 
tion on magnetic substrates [49], and correctly assessing 
formation energies m- While the nonlocal functional 
VV09 considers spin in its own way [SU [52], there have 
so far only been pragmatic approaches for vdW-DF— 
ignoring the effect of spin on the nonlocal correlation al¬ 
together [53] or estimating the effect [541156] using the 
semi-local correlation of PBE [28]. 

In this letter, we formulate a proper extension of vdW- 
DE to spin-polarized systems, termed svdW-DE, follow¬ 
ing the design-logic of the original functional. We apply 
svdW-DE to study the nonlocal-correlation driven ad¬ 
sorption of H2 and GO2 in MOE74 and find that spin 
plays a significant role, providing a detailed analysis of 
spin signatures in such vdW bonding. Beyond MOEs, 
we envision that svdW-DE will lead to wider materials- 
theory progress in a stimulating role like that of LSDA, 
i.e., LDA’s spin extension m- LSDA was introduced 
to describe bulk-cohesive and molecular-binding energies 
[5Q[ [57I [58] but also led DET to important successes in 
the study of magnetism [59]. The svdW-DE formulation 
enables a robust exploration of systems where spin and 
nonlocal correlations are both important and it makes 
vdW-DE a general purpose method ED]. 

To design svdW-DE as the natural extension of vdW- 
DE to spin-polarized systems, we revisit the derivation 
of its nonlocal correlation energy functional. The start¬ 
ing point is the adiabatic-connection formula (AGE) ex¬ 
pressed in terms of a scalar dielectric function e that re¬ 
flects a formal average of the coupling-constant integra¬ 
tion over the screened density-response function [231 [sa 
ED Eg. This provides a split-up of the total exchange- 
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correlation energy into a nonlocal and semi-local piece 
defined by [22[ [23l [25l [30] in terms 
of the Coulomb Green’s function G and an integral over 
the imaginary frequency u: 

d?/ 

‘ [n]= Tr [ In (Ve ■ VG) - In e] (1) 

Jo 27r 

pOC 1 

^'"c[n] = J^ ^Trlne-E^^u. ( 2 ) 

The general-geometry vdW-DF versions musiiisi ex¬ 
pand the nonlocal correlation energy 0 in terms of a 
semi-local response function S = \ne that is parame¬ 
terized via the choice of internal semi-local (GGA-type) 
functional E^^^ ([^. 

To obtain a computationally tractable approximation 
for E^^ [22] |23l [25], vdW-DF relies on a plasmon-pole 
approximation of S defined in plane-wave representation 
as with 


/ 
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Here, n(r) is the total electron density and cjg(r) is 
the effective local plasmon dispersion, parameterized by 
an effective response parameter in the form of an in¬ 
verse length scale go(r) = qo[n] = go(^(r), Vn(r)) ; m 
and e are the electronic mass and charge. The link 
between the energy-per-particle of the internal func¬ 
tional (r) and ci;g(r) [22] |23[ [30] follows from com¬ 
bining Eqs. 0 and 0 together with a plasmon disper¬ 
sion co’g(r) = g^/2h(g/go(i*)) with a Gaussian shape of 
h{x) = 1 — exp(—where 7 is an arbitrary constant 
set to 47 r/ 9 . Expanding Eq. 0 to second order in S', one 
arrives at the well-known six-dimensional integral over a 
universal kernel 4>o(a,6), 


Ef = i y drdr' n(r) $o(?o(r)|r-r'|,go(r')|r-r'|) n(r') , 

(4) 

which defines the approximation for E^^ [23]. The to¬ 
tal exchange-correlation energy also consists of 

the semi-local functional E^^^ (0. This is in practice ap¬ 
proximated as E^^^ ^ E^^ = based on a 

number of criteria 1211 ia Eli ED isa [S 3 and differing 
from the internal functional E^^^ to varying degrees. 

The extension of the semi-local part E^^ to spin- 
polarized systems is straightforward. It is given by the 
exact spin scaling of exchange [63], i.e. = 

Ex[2n^]/2 + Ex[2ni]/2, and the well-established spin- 
dependence of the local correlation [27]. Here and 
denote the spin-density components. Grucially, by apply¬ 
ing the very same criteria, we obtain a fully consistent 
extension of E^^ for the spin case. 

The spin-scaling of exchange results in a spin- 
dependent semi-local response S in Eq. 0 , with spin 
entering exclusively in the denominator through uJq{T). 


The numerator is given by the /-sum rule, specified as the 
classical plasmon frequency which depends only on the 
total electron density n(r). The formulation of svdW-DE 
can therefore be based on an universal-kernel evaluation 
using the exact same function 4>o(a,6) as in vdW-DE. 
Nevertheless, the form of the effective response parame¬ 
ter go—which acts as scaling parameter in the arguments 
of 4>o—must be adjusted, qo[n] go[^t 5 %]’ reflect the 
explicit spin dependence of the plasmon dispersion. 

Motivating our procedure for extending the original 
vdW-DE formulations to spin-polarized system is the 
interpretation of the vdW-DE nonlocal correlation en¬ 
ergy as a formal summation of zero-point energy shifts 
[3Ql [ 6 ^ [65] . The vdW-DE framework starts with a de¬ 
scription of the semilocal exchange-correlation holes cor¬ 
responding to the internal functional E^^^ [22] [25] [29] . 
using a plasmon model to characterize the associated re¬ 
sponse. The vdW-DE nonlocal correlation energy 0 
is a rigorous summation of the plasmon-pole shifts that 
result when such holes couple electrodynamically m- 
Spin clearly affects the GGA-type internal hole and our 
svdW-DE formalism represents a proper implementation 
of how such semilocal spin effects impact the summation 
of zero-point energy shifts in Eq. 0 . 

To establish the updated form of go['^t’^t] 
structive to first revisit how go [n] is specified in the spin- 
neutral case, where it is given as scaling of the Eermi 
wave vector /ci 7 (r) = ( 37 r^n)^/^ as follows 

= lda/ = qoc[n] + qoxM , (5) 

4'7F 

gocM = -^7°^, (6) 

qox[n] = -(l - 04“"^ • (7) 

Here and the exchange gradient cor¬ 

rections are expressed in terms of the scaled gradient 
5 = |Vn|/2/ci?n. These relations 0-0 can be directly 
adapted to the spin case. The correlation part goc[^t 5 %] 
is specified by the spin-dependent PW92 EDA correla¬ 
tion energy-per-particle [27] . Eor the exchange part 

qox[^\^'^\\i the spin scaling relation [63] gives the follow¬ 
ing form 


?T/T 77/1 

p-^goa;[2%] + ^^_^^^ gOa;[2n4.] . (8) 


These equations fully specify the nonlocal correlation en¬ 
ergy of svdW-DE. We make svdW-DE self-consistent, im¬ 
plemented in Quantum Espresso [ 66 ], by computing 
the corresponding exchange-correlation potential [24] . 
Eurther details on the implementation and calculations 
are provided in the Supplemental Material. 

We test svdW-DE on three cases of increasing com¬ 
plexity; results are summarized here and details are in 
the Supplemental Material. Eor our test cases we use 
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the implementations svdW-DFl [23] and svdW-DF2 [25| 
(which are better suited for small molecules) as well as 
svdW-DF-cx [26] (which is better suited for larger, ex¬ 
tended systems). We start with the Li-dimer in its triplet 
state —an ideal test case that critically balances vdW 
and spin effects. We find a dissociation energy of 53 meV 
for svdW-DFl and 70 meV for svdW-DF-cx, the for¬ 
mer in good agreement with the experimental value of 
41 meV [72]; VVIO and PBE find 77 meV. A second case 
is given by atomization energies for molecules from the 
G1 set ESI, where spin enters through magnetic molecu¬ 
lar ground states and the isolated atoms. We find a mean 
absolute percentage error of 4.59% and 7.75% for svdW- 
DFl and svdW-DF-cx; VVIO and PBE find 5.14% and 
7.11%, respectively. Eor a third, extended-system test, 
we study the weak-chemisorption of graphene on Ni(lll) 
[74H76], finding a binding separation for svdW-DE-cx of 
2.12 A, in excellent agreement with experiment (2.11 ± 
0.07 A [74]). In contrast, svdW-DEl finds 3.76 A, VVIO 
finds 3.37 A, and PBE essentially does not bind—unlike 
svdW-DE-cx they all miss a significant chemical compo¬ 
nent to the binding. 

Table [I] summarizes the main point of this letter: that 
svdW-DE provides insight on the nature of nonlocal spin 
effects in the adsorption of H 2 and CO 2 in the linear 
magnets Mn-MOE74, Ee-MOE74, Co-MOE74, and Ni- 
MOE74 [77]. The table reports raw svdW-DE binding 
energies AE, as well as values AEzpe corrected for zero- 
point vibrations of the adsorbates and binding enthalpies 
at room temperature Ai7298- We note that the adhesion 
comes entirely from —without nonlocal correlations, 
CO 2 would not bind at all and H 2 would only bind with 
a binding energy of ^5 meV. 

Overall, we find very good agreement with experiment, 
which we partly attribute to the “ex” version of svdW- 
DE; agreement with other vdW-DE calculations is also 
good [ 5 ]. In all cases we find that the inclusion of spin 
has an important effect on the binding. Note that the 
case of Mn is somewhat artificial and the “no spin” num- 


TABLE I. Binding energies [meV] of small molecules in the 
system A4-MOF74-hA with M = Mn, Fe, Co, and Ni and A 
= H 2 and CO 2 . In number triplets the first number refers to 
the bare binding energy AE, the second one includes the zero- 
point correction AEzpe, and the third refers to the binding 
enthalpy at room temperature Ai7298- 


M 

A 

Exp. 

no spin 

spin 

Mn 

H 2 

91 EZ] 

-473/-524/-524 

-133/-117/-117 

Fe 

H 2 

104 EH] 

-134/-137/-137 

-122/-124/-124 

Co 

H 2 

111 Eli 

-177/-181/-181 

-111/-117/-117 

Ni 

H 2 

134 Ell 

-134/-137/-137 

-131/-133/-133 

Mn 

CO 2 

331 ESI 

-528/-551/-550 

-337/-345/-344 

Fe 

CO 2 

352 EH] 

-344/-353/-351 

-315/-323/-321 

Co 

CO 2 

383 [70] 

-377/-387/-385 

-350/-359/-357 

Ni 

CO 2 

394 [n] 

-300/-311/-309 

-377/-390/-387 



spin up spin down 


FIG. 1. (Upper panels) Induced charge density upon CO 2 
adsorption in Ni-MOF74. (Lower panels) Induced charge den¬ 
sity split into its up and down contribution. Blue (yellow) 
areas show charge depletion (accumulation). Iso levels are 
0.001 e/Bohr^. See Fig. I in the Supplemental Material for 
the structure of the MOF. 


bers seem inflated—this is a result of the fact that Mn- 
MOF74 itself requires spin for a proper description of its 
structure. 

Of particular interest is the binding of GO 2 in Ni- 
MOF74, where spin effects play a tantalizing and unex¬ 
pected role. The GO 2 molecule is diamagnetic and should 
experience a slight repulsion and weaker binding in the 
presence of the magnetic dipole of Ni—similarly to what 
is observed in all the other cases in Table [I| However, on 
the contrary, when spin effects are included the binding 
increases and the molecule experiences a stronger attrac¬ 
tion, which warrants further investigation. In the upper 
panels of Fig. we plot the induced charge density, i.e., 
the charge density redistribution due to the formation of 
the bond. It is clearly visible that in the spin case more 
charge is pulled in-between the GO 2 and the metal site, 
resulting in the stronger binding. In the spin case, we 
can split this induced charge density further into spin-up 
and spin-down contributions, as shown in the lower pan¬ 
els of Fig. Here we see the true spin effect: much more 
spin-down density is being pulled into the bond, com¬ 
pared to spin-up density. This peculiar behavior can be 
understood by analyzing the projected density of states 
in Fig. In particular, from the middle and bottom 
panel we see that at — 5 eV the O p states show similar 
peaks in the spin up and down channels. However, the 
projected Ni d states at that point have a large spin down 
density while the spin up density is much smaller. Thus, 
the O p states hybridize with the Ni spin-down d states. 
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E-Ep [eV] 


FIG. 2. Up (black) and down (red) total density of states 
(top) and projected density of states on the Ni d states (mid¬ 
dle) and O p states (bottom) of the Ni-M0F74+C02 system. 


while the hybridization with the spin-up states is negligi¬ 
ble (see Fig. in the Supplemental Material for plots of 
the corresponding orbitals). The interaction of the O p 
states with the down-spin Ni d states is therefore respon¬ 
sible for the increased and counter-intuitive strength of 
the bond. 

As mentioned above, is responsible for the entire 
binding of the small molecules. As such, it is at least 
indirectly responsible for all the effects we see in our 
figures and tables. To examine the spin effect of 
explicitly, we also calculate the difference between the 
binding-induced density of a full svdW-DF calculation 
and the same calculation without the E"^^ term (Fig. 
Supplemental Material). Although overall smaller in 
magnitude—as expected, since the semi-local part E^^ 
also contributes to the induced density—we find the same 
behavior as in Fig.[2 more down density is being pulled 
into the bond, strengthening the binding. 

It is also revealing to partition the charge and magnetic 
moment of the system, as detailed in Table [III] in the Sup¬ 
plemental Material. Partitioning schemes are not unique, 
but one can still gain qualitative information. Before ad¬ 
sorption, the CO 2 has no magnetic moment and all six Ni 
atoms in the unit cell are equivalent. However, once ad¬ 
sorption occurs, the up and down charge inside the CO 2 
rearranges differently and gives rise to a small but observ¬ 
able magnetic moment. At the same time, the adsorp¬ 
tion process leaves the total charge on the CO 2 molecule 
unchanged, i.e. there is no net charge transfer. During 
this adsorption process, the nearby Ni loses 0.03 e—our 
analysis shows that it loses mostly down density—and 
thus its magnetic moment increases, giving rise to a weak 
nonlocal-correlation induced magnetic interaction. 

Finally, in terms of absolute numbers, spin effects 
in the Fe-M0F74+C02 and C 0 -MOF 74 +CO 2 systems 


seem to worsen agreement with experiment—but, at the 
same time, they actually resolve a more pressing issue. 
Experimentally, the CO 2 binding strength should follow 
the order Mn < Fe < Co < Ni. In calculations without 
spin (not considering the artificial case of Mn), the or¬ 
der is reversed. However, after including spin effects the 
correct order is restored. 

In summary, we have developed a consistent spin- 
polarized version of the nonlocal exchange-correlation 
functional vdW-DF, which we find to now become an 
all-purpose functional. We then apply this framework 
to study small-molecule adsorption in MOF74 with mag¬ 
netic open metal sites and find that including nonlocal 
spin effects can significantly influence the binding of the 
adsorbates. In the case of Ni-M0F-74+C02 we find a 
counter-intuitive increase in binding due to nonlocal spin 
effects, which we explain by a coincidental interaction of 
the Ni d and O p states. The additional degree of freedom 
from such unexpected magnetic interactions can be used 
to tailor the specificity of MOFs in novel gas storage, 
sequestration, and sensing applications. 

Work in the US has been supported by NSF Grant 
No. DMR-1145968 and DOE Grant No. DE-EG02- 
08ER46491. Work in Sweden has been supported by 
grants from the Swedish Research Council (VR), the 
Swedish Eoundation for Strategic Research (SSE), the 
Chalmers e-Science Centre, and the Chalmers Materials 
Area of Advance. 
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— Supplemental Material — 


I. COMPUTATIONAL DETAILS 

We have implemented svdW-DF in PWSCF 5.1, which 
is part of the Quantum Espresso package [ 66 ]. We 
used ultrasoft pseudopotentials and a plane wave cutoff 
for the wave functions of 40 Ryd together with the ex [26] 
implementation of svdW-DF. Hubbard U corrections are 
applied to the d electrons of Mn, Fe, Co, and Ni, us¬ 
ing the well-established values of 4.0, 4.0, 3.3, and 6.4 eV 
from Ref. [78] , where the U values for each metal A4 were 
determined so that they reproduce the experimental oxi¬ 
dation energy of the metal monoxide to AI 2 O 3 . Note that 
these 1/ values have been shown to give highly accurate 
results for small-molecule adsorption in MOFs U SO]. 
Due to the large unitcell size of the MOF, we only sam¬ 
pled the F-point. In all our calculations, we optimized 
the entire structure until the forces on all atoms reached 
less than 1 meV/A. 

The starting points of the calculations are the experi¬ 
mental rhombohedral structures of Mn-, Fe-, Co-, and 
Ni-MOF74 with 54 atoms in the primitive cell. Ta¬ 
ble [TI] shows the experimental lattice constants, listed 
in the more convenient hexagonal lattice [671 [79]. For 
the magnetic ordering of the metal ions we adopted the 
ferromagnetic arrangement, which has been found to be 
the ground-state configuration m Figure provides 
a graphical representation of the isostructural Mn-, Fe-, 
Co-, and Ni-MOF74. 

Vibrational properties were calculated using the finite 
difference method with atomic displacements of 0.015 A 
in each direction. The corresponding Hessian matrix was 
symmetrized and the acoustic sum rule was enforced. 
Only atoms of the adsorbed molecule were allowed to 
move—an approximation which we found to result in er¬ 
rors on the order of merely 1 cm“^. The zero-point en¬ 
ergy and the thermal correction to the binding energy 
was calculated as 


AiJzPE + AiJr = 9 E E -?—Tibrn —T ’ 

2 ^ ^ exp{vi/kBT) - 1 

where is Boltzmann’s constant and Vi are the vibra¬ 
tional frequencies of the system in eV. The temperature 
was set to T = 298 K. Contributions from low frequency 
vibrations were not included. 


TABLE II. Experimental lattice constants used [A] for our 
A4-MOE74 systems with A4 = Mn, Ee, Co, and Ni. 


M 

Ref. 

a = b 

c 

Mn 

m 

26.230 

7.035 

Ee 

m 

26.098 

6.851 

Co 

m 

25.948 

6.838 

Ni 

m 

25.719 

6.741 



EIG. 3. MOE74 structure with a CO 2 molecule adsorbed at 
one of the open metal sites. The hexagonal channel structure 
with its six equivalent open metal sites per unit cell is clearly 
visible. Carbon atoms are depicted as grey spheres, while oxy¬ 
gen, hydrogen, and metal atoms are shown in red, white, and 
blue. The rectangle indicates the portion of MOE74 shown in 
the main manuscript and Eigs.^and[^ 
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II. ADDITIONAL TABLES AND FIGURES 



TABLE III. Bader charge q (see Refs. [80|[HI]) [in units of e, 
showing electron loss relative to the neutral unit] and mag¬ 
netic moment fi [in units of jis] during the CO 2 adsorption 
process in Ni-MOF74. * indicates the Ni atom participat¬ 
ing in the bond. Only one representative of the remaining 5 
equivalent Ni is given. 



before adsorption 

after adsorption 

Unit 

q 


q 


CO 2 

0.00 

0.00 

0.00 

0.01 

Ni* 

1.30 

1.57 

1.33 

1.60 

Ni 

1.30 

1.57 

1.30 

1.57 


O 

2 



spin up 


spin down 


FIG. 4. Contour plots of |'0|^ x sign('0) evaluated (at 
the F-point) for the spin-up and spin-down Ni-MOF74 bands 
around — 5 eV with (top row) and without (bottom row) the 
adsorbed CO 2 molecule. Iso levels are 0.0002 e/Bohr^. In all 
cases, the d^ 2 -\\k.e orbital on the Ni is visible. In the spin-up 
channel, the adsorption of CO 2 leaves this orbital almost un¬ 
changed. In contrast, in the spin-down channel, the dz 2 -Y\]^e 
orbital rotates towards the oxygen p orbital to form a new, 
hybridized orbital. 



spin up 


spin down 


III. SELF-CONSISTENT DERIVATIVES 


As for the evaluation of goc[^t,^t] derivatives, 

they are completely defined through the choice of cor¬ 
relation inside and are given elsewhere [27] • For 

follows that 


dqox[n^,nx\ 

dnt 


^4 


+ 


(nt + n4,)2 
2nt dqox[n 


( 90 a:[2%] - qox[‘^rn\) 


dn 

dqox[nt,n^] _ 


n=2n>|' 


(lOa) 


dn. 


+ 


(n-f + n4,)2 

2^4, dqox[n 


{qox[‘^ni\ - qox[2nf]) 


n=2ni 


and 


dn 

dqox[n^,ni,] ^ n-|- dqox[2n^] 

d\Vn^\ Uf + ni (i|Vn-|-| 

dqox[n^,ni] ^ ni dqox[2n^] 

d\Vni\ Uf + n^, d|Vn^| 


(10b) 

(lla) 

(llb) 


FIG. 5. Induced charge density originating from the 
term upon CO 2 adsorption Ni-MOF74, split into its up and 
down contribution. Blue (yellow) areas show charge depletion 
(accumulation). Iso levels are 0.0001 e/Bohr^. 
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IV. TESTING SVDW-DF 

Testing and benchmarking has been performed on three model systems, i.e. the Li dimer in its triplet state 
atomization energies of a set of small molecules, and graphene on the Ni(lll) surface. We used our implementation 
in PWSCF [66] and calculated binding curves and energies using svdW-DFl [23|, svdW-DF2 [25], svdW-DF-cx [26] . 
VVIO [38|[82], and PBE [28]; these studies were performed using pseudopotentials from the Rutgers database with the 
recommended cutoffs [83]. Note that svdW-DFl and svdW-DF2 are better suited for small molecules, while svdW- 
DF-cx is better suited for larger, extended systems. We compare our PWsCF calculations to high-level quantum 
chemistry (QC), diffusion Monte Carlo (DMC), and experiment where appropriate. Please see the footnotes [51] and 
[53] in the main text concerning the neglect of spin-polarization effects in the nonlocal part of VVIO. 


A. Li Dimer in the Triplet State 


TABLE IV. Dissociation energy De [meV] and fundamental frequency Ue [meV] for Li 2 in the triplet state ^E. The experimental 
value is =41 meV [72] [84]. Early papers using quantum chemistry methods report a considerable spread in De values 

[85h87j . but the most recent high-level quantum chemistry (QC) results show in essence perfect agreement with experiment [48] . 



svdW-DFl [13] 

svdW-DF2 US] 

svdW-DF-cx [26] 

VVIO 131] 

PBE Hi] 

QC 131] 

De 

53 

55 

70 

77 

77 

41 

CUe 

8 

8 

8 

8 

8 

8 







B. Atomization Energy of Small Molecules 


TABLE V. Atomization energies Eat [eV] of small molecules from the G1 set [73]. Spin expectation values {S{S + 1)) different 
from zero indicate a magnetic ground state. Diffu sion Monte Carlo (DMC) values are taken from |88j : experimental numbers 
from [89]. An error analysis follows in Table VI Spin enters through the isolated atoms and/or magnetic molecular ground 
states. 



(5(5+1)) 

svdW-DFl ES] 

svdW-DF2 EH 

svdW-DF-cx EH 

vvio EHl 

PBE EH 

DMC EH 

Exp. eh 

CHa 

3/4 

13.421 

13.363 

13.639 

13.373 

13.507 

12.615 

12.545 

CH4 

0 

18.074 

17.937 

18.502 

17.963 

18.212 

17.129 

17.020 

NH 

2 

3.763 

3.810 

3.758 

3.795 

3.778 

3.393 

3.426 

NH2 

3/4 

7.942 

7.934 

8.100 

7.957 

8.042 

7.339 

7.372 

NHs 

0 

12.648 

12.516 

13.088 

12.620 

12.880 

11.991 

11.999 

OH 

3/4 

4.236 

4.098 

4.557 

4.195 

4.408 

4.390 

4.397 

H2O 

0 

9.660 

9.449 

10.194 

9.612 

9.953 

9.514 

9.512 

SiHs 

3/4 

9.887 

9.863 

10.015 

9.594 

9.668 

9.328 

9.280 

SiH 4 

0 

13.921 

13.941 

14.086 

13.528 

13.595 

13.261 

13.122 

PH2 

3/4 

6.716 

6.694 

6.899 

6.572 

6.656 

6.232 

6.275 

H2S 

0 

7.901 

7.786 

8.190 

7.871 

8.045 

7.463 

7.506 

HCl 

0 

4.575 

4.457 

4.785 

4.575 

4.697 

4.486 

4.432 

LiE 

0 

5.885 

5.990 

5.984 

5.897 

6.209 

6.292 

5.984 

C2H4 

0 

24.272 

24.042 

25.020 

24.439 

24.819 

23.136 

23.065 

CO 

0 

10.940 

10.732 

11.489 

11.179 

11.501 

10.981 

11.110 

N2 

0 

9.590 

9.494 

10.017 

9.972 

10.164 

9.587 

9.761 

NO 

3/4 

6.373 

6.133 

7.043 

6.725 

7.079 

6.198 

6.507 

O2 

2 

5.171 

4.805 

6.017 

5.434 

5.925 

4.844 

5.115 

CO2 

0 

16.526 

16.028 

17.768 

17.103 

17.796 

16.458 

16.562 

Na 2 

0 

0.794 

0.767 

0.740 

0.671 

0.707 

0.750 

0.729 

S2 

2 

4.820 

4.650 

5.197 

5.161 

5.334 

4.264 

4.365 

NaCl 

0 

4.085 

3.993 

4.190 

4.129 

4.157 

4.284 

4.219 

SiO 

0 

8.045 

7.971 

8.508 

8.229 

8.429 

8.096 

8.239 

cs 

0 

7.504 

7.374 

7.861 

7.801 

7.980 

7.172 

7.328 

so 

2 

5.589 

5.357 

6.156 

5.866 

6.178 

5.100 

5.351 


TABLE VI. Error analysis for the datasets given in Table [vj Given are the mean (signed) percentage error (MPE), mean 
absolute percentage error (MAPE), and mean absolute deviation (MAD) values when comparing with experiment. Note that 
svdW-DE-cx, due to its design, is not necessarily expected to perform well for small molecules, but still gives good results, 
comparable to PBE. 



svdW-DFl EH 

svdW-DF2 EH 

svdW-DE-cx [ 2 ^ 

vvio EH 

PBE EH 

DMC eh 

MPE [%] 

-3.28 

-1.52 

-7.69 

-3.84 

-6.75 

0.54 

MAPE [%] 

4.59 

4.50 

7.75 

5.14 

7.11 

1.65 

MAD [eV] 

0.37 

0.38 

0.67 

0.40 

0.59 

0.11 
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C. Graphene on Ni(lll) 



FIG. 6. Unit cell used for the adsorption of graphene on 
the Ni(lll) surface. Ni was modeled as a slab consisting of 
6 layers of atoms, the bottom three of which were hxed to 
their bulk positions. Blue and grey spheres denote Ni and C 
atoms, respectively. 



Distance [A] 


FIG. 7. Binding energy curve of graphene on the Ni(lll) 
surface. The dashed vertical line denotes the experimental 
value for the distance between graphene and the Ni(lll) sur¬ 
face [74] . 
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